Small mass- and trap-imbalanced two-component Fermi systems 
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Motivated by the prospect of optical lattice experiments with two-component Fermi gases con- 
sisting of different atomic species such as Li and K, we calculate the energies for N fermions under 
harmonic confinement as a function of the mass- and trap-imbalance, i.e., as a function of the ratio 
between the masses and frequencies of species one and two, using microscopic approaches. Our en- 
ergies for N — 2 through 6 can be used to determine the energetically most favorable configuration 
for a given number of atoms per species of a deep lattice in which each lattice site is approximately 
harmonic and in which tunneling between neighboring sites can be neglected. Furthermore, our 
energies determine one of the input parameters, namely the onsite interaction strength, of the corre- 
sponding lattice Hamiltonian. We also determine and interpret the excitation gap for unequal-mass 
systems with up to N = 13 atoms for equal oscillator lengths. 

PACS numbers: 



I. INTRODUCTION 

Cold-atom experiments have reached an impressive 
level of sophistication over the past decade. About ten 
different atomic species have been Bose condensed and, 
although experimentally more challenging, an increasing 
number of fermionic species have been cooled to quan- 
tum degeneracy, including 3 He [lj], 6 Li 0, 40 K Q and 
two Yb isotopes ( m Yb and 173 Yb) 0, [J To date, ex- 
periments on fermionic atoms have focused on studying 
Bose-Fermi mixtures 0, 0, H, HI, one-component Fermi 
systems with p-wave interactions [lOj and equal-mass 
two-component Fermi systems with interspecies s-wave 
interactions [HE2E1- 

Presently, the simultaneous trapping and cooling of 
two different fermionic species is being actively pursued 
by a number of laboratories [3, El, EH , adding a new 
degree of freedom, i.e., the mass ratio between the two 
atomic species. Unequal-mass two-component Fermi sys- 
tems are expected to behav e q uite differently than the 
equal-mass counterpart [13, El EIS HH |M El, [H, 
HE HI]. From the few-body perspective, the existence 
of weakly-bound trimers for sufficiently large mass ratios 
consisting of two heavy fermions and one light fermion is 
intriguing [Tt], El, EiJ HI] ■ Whether the existence of these 
bound trimer states allows, e.g., for the formation of a 
gas consisting of trimers with sufficiently long lifetime 
has been discussed [13] ■ On the other hand, adopting a 
many-body perspective [U, [12 [H, Hi| , the ground state 
phase diagram of mass- and population-imbalanced two- 
component Fermi systems has been predicted to show 
quantum and topological phase transitions which are 
not present in the phase diagram of population-balanced 
equal-mass two-component Fermi systems. 

The increasing interest of not only the atomic physics 
community but also the nuclear physics, molecular 
physics, condensed matter physics and quantum infor- 
mation science communities in cold atom systems can be 
attributed to two major achievements. First, the atom- 
atom scattering length can be adjusted experimentally to 



essentially any value, including vanishingly small and in- 
finitely large positive or negative values, by applying an 
external field in the vicinity of a so-called Fano-Feshbach 
resonance [H, [29|, [H, [3l| . Second, cold atomic gases can 
be loaded into an optical lattice [12 [H, [H, HH , allowing, 
e.g., for the study of the Mott-insulator transition [32j], a 
topic historically primarily considered by condensed mat- 
ter physicists. Furthermore, cold atom systems loaded 
into optical lattices may ultimately be used as a quan- 
tum simulator [H, H3, HI] . 

To date, most microscopic studies of equal-mass sys- 
tems have assume d eq ual trapping p otentials of the two 
species [U S3, ED, 52 0, 51 El El- However, the lat- 
tice potential felt by the two different hyperfine states 
may be different even for equal-mass systems, leading 
to trap-imbalanced systems (47|. For unequal-mass sys- 
tems such as a 6 Li- 40 K mixture (for which the mass ra- 
tio k is approximately 6.7), the trapping potentials felt 
by the two species are, in general, different, owing to 
the mass difference and the species-dependent proper- 
ties of the hyperfine states. The trapping potentials felt 
by the two species may be tuned to some degree exper- 
imentally [2J, 53- Motivated by these considerations, 
the present paper explores the rich behavior of trap- and 
mass-imbalanced systems. These systems share some 
similarities with population-imbalanced systems [H, 5t| , 
which have received considerable attention recently. 

Assuming a deep lattice with neglegibly small tunnel- 
ing between neighboring lattice sites, this paper deter- 
mines the ground state properties of small s-wave inter- 
acting two-component Fermi systems trapped by spheri- 
cally symmetric harmonic species-specific potentials with 
trapping frequencies u\ and u>2, respectively. Through- 
out, we adopt a microscopic many-body framework. Our 
main results are: (i) For small but negative s-wave scat- 
tering lengths, we determine a compact expression for 
the ground state energy of small systems with unequal 
masses and trapping frequencies perturbatively. (ii) In 
the strongly-interacting unitary regime, our numerical 
energies determine the phase diagram of optical lattice 
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systems in the no-tunneling regime for a large range of 
mass ratios and trapping frequencies and the on-site in- 
teraction strengths that parametrize the corresponding 
lattice Hamiltonian; furthermore, they provide insights 
into the behavior of the excitation gap. (iii) We show 
explicitly that the behavior of trap-imbalanced systems 
with small and positive s-wave scattering lengths is — just 
as that of trap-balanced systems — to leading order gov- 
erned by the dimer-dimer scattering length. 

Section HT1 introduces the Hamiltonian and the numeri- 
cal techniques employed to solve the corresponding time- 
independent Schrodinger equation. Section IIIII contains 
our results for small negative, infinitely large and small 
positive s-wave scattering lengths. Finally, Sec. |IV] con- 
cludes. 



II. THEORETICAL BACKGROUND 

A. Hamiltonian 

The adopted model Hamiltonian H for a two- 
component Fermi system with N\ mass mi and N2 mass 
mi atoms under spherically harmonic confinement reads 

i— 1 x 7 

j'=l v z ~ ' i=l i'=l 

where fi and ry denote the position vectors of the ith 
atom of species 1 and the i'th atom of species 2, respec- 
tively, and 0J1 and u>2 the angular trapping frequencies 
felt by the atoms of species 1 and 2, respectively. The 
interaction potential V depends on the interparticle dis- 
tance vector fn> , fai = fi — t> , and is characterized by 
the s-wave scattering length a s . Throughout, like atoms 
are assumed to be non-interacting, which is well justified 
for most experimentally relevant systems. 

Our perturbative, small \a s \ analysis (see Sec. JIII A]) 
considers a zero-range S- function potential V$(f) [50(, 

V 5 (f) = 2 -^8{r), (2) 
A 1 

where fi denotes the reduced mass, /1 = mims/fmi + 
7712). In our numerical calculations (see Sees. llTTBl and 
III CP we employ, as in our previous calculations [4J, |45| . 
46| , a shape-dependent spherically-symmetric square well 
potential V sw (r) with range R and depth Vq (Vq > 0), 

V ™( r ) = { forr>i? ' ® 

where r = \r\. For a fixed Rq, Vq is adjusted so that the 
interspecies s-wave scattering length a s takes on the de- 
sired value. Section fill Bl considers the so-called unitary 



regime, where Vq is adjusted so that the two-body poten- 
tial supports a zero-energy s-wave bound state, implying 
a diverging s-wave scattering length a s , i.e., \/a s = 0, 
but no deeply-lying bound states. Section IIII CI in 
contrast, considers the regime where Vq is adjusted so 
that the free-space dimer supports one deep-lying s-wave 
bound state (whose binding energy depends on the de- 
tails of the two-body potential) , implying a small positive 
s-wave scattering length. The range Rq of V sw is taken 
to be small compared to the oscillator lengths aho.i, 



aho,i = \/h/(m,iU)i), (4) 

where i = 1 or 2. Most calculations reported below use 
Rq = 0.01a^ Oj i. To estimate how well the resulting prop- 
erties agree with those for zero-range interactions, we 
analyze the dependence of the observables on the range 
Rq in detail for a few selected cases. 

Section HITl presents our results for three different scat- 
tering length regimes, i.e., for weakly-attractive Fermi 
gases (\a s \ small and a s < 0), for strongly-interacting 
Fermi gases (l/|a s | = 0) and for weakly-repulsive Fermi 
gases (a s small and a s > 0). In all three regimes, we de- 
termine the energies of small trapped systems with either 
Ni = N2 or \N\ — JV2I = 1. In addition to changing the 
number of particles and the scattering length a s , we vary 
the mass ratio k, 

« = m2/mi, (5) 

and the ratio LO2/0J1 between the two trapping frequen- 
cies. Mass ratios ranging from k — 1 to 8 are consid- 
ered (for unequal-mass systems species 2 has the heavier 
mass). For k > 8.6, three-body bound states have been 
predicted to exist for systems that consist of two heavy 
fcrmions and one light fcrmion and that interact through 
zero-range potentials pjj [2(| • While studying the impli- 
cations of these three-body states for many-body systems 
is interesting (see, e.g., Ref. |13|), this topic is beyond the 
scope of the present paper. 



B. Numerical techniques 

To solve the time-independent Schrodinger equation 
for the Hamiltonian given in Eq. ([1]), we employ two dif- 
ferent numerical techniques. For Ni = N2 = 1, we first 
build and then diagonalize the Hamiltonian matrix while 
we resort to the fixed-node diffusion quantum Monte 
Carlo (FN-DMC) technique p5ll . |52| for larger systems. 

We first discuss the diagonalization approach employed 
to solve the Schrodinger equation for the Hamiltonian 
given in Eq. (TT]) with N± = N 2 = 1 and V = V sw ; it 
follows Ref. [53| , with the main difference that we use 
a finite-range square-well potential while Ref. [53[ uses 
the Fermi-Huang pseudo-potential [54| . We rewrite our 
two-body Hamiltonian H t b in terms of a center-of-mass 
Hamiltonian H cm , a relative Hamiltonian H re i and a cou- 
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pling term V coup (R,r) [Hj], 

H tb = H cm + H rel + V coup (R, r) , (6) 

where 

Hem = + \Mu? cm R\ (7) 

Hrel = ^-^l+\lM4eir\ (8) 

and 

V COU p(R, r) = ^lou P R ■ r. (9) 

Here, R and r (r = r\ — r*2) denote the center-of-mass 
and relative vectors, respectively, and M denotes the to- 
tal mass of the two-body system, M — m\ + m^- The 
frequencies Lo cm , uj re i and uj coup are defined as 



w™ = V [miwl + m 2 w|)/M, (10) 



U re \ = sj {miUj\ + miJ$)/M (11) 

and 



Wcoup = y \ui ~ (12) 

For equal trapping frequencies, vanishes and 

both w cm and oj re i reduce to wi (which equals lu 2 ). In 
this case, the center-of-mass and relative motions decou- 
ple, and the total wave function ^(R,r) can be written 
as a product of a independent function <$>nlm l an d a 
r-dependent function (f> n i mi . The &nlm l and <fi n lmi are 
solutions to the Schrodinger equations for H cm and H re i , 
respectively, and the subscripts NLM L and nlmi denote 
the principal, angular momentum and projection quan- 
tum numbers of the center-of-mass and relative systems, 
respectively. The &nlm l are the harmonic oscillator 
wave functions of a mass M particle with eigenenergies 

E NL = (2N + L+^J?tiJ cm , (13) 

where N = 0, 1, • • • , L = 0, 1, 2, • • ■ and M L = —L, -L + 
1, • • • , L. For the spherically-symmetric square well po- 
tential V sw , the angular part of the relative wave function 
4>nimi is given by the spherical harmonic Yi mi while the 
radial part R n i can be written in terms of the confluent 
hypergeometric function M for r < Rp and the Kummer 
function U for r > Ro (see, e.g., Ref. 55]). Equating the 
log-derivative of the inner and outer radial wave func- 
tions at r = i?o results in a compact expression for the 
eigenequation, from which we obtain the eigenenergies 
E„i of H re i using standard root-finding techniques. The 



radial wave functions R n i are then readily obtained by 
enforcing continuity at r — Rq. We normalize the i?„;(r) 
numerically. 

To determine the eigenenergies of the two-particle 
Hamiltonian H t b with non-zero V coup (R, f), we expand 
the full wave function ^(R, r) in terms of the complete set 
{®NLMi,(R)<f>nlmi(r)}- Recognizing that Htb commutes 
with the z-component of the total angular momentum 
operator (i.e., that Ml + mi is conserved), we restrict 
the allowed Ml and mi combinations to Ml + mi = 0. 
Since the &nlm l and <fi n lmi are solutions of H cm and 
H re i , respectively, H cm and H re i are diagonal in this rep- 
resentation. To evaluate the matrix elements involving 
V coup , we rewrite the dot product R ■ r in terms of R, 
r, and the spherical harmonics Yl=\,m l an d Yt=i,m t as- 
sociated with the center-of-mass and relative degrees of 
freedom, respectively. The angular integrals then read- 
ily reduce to Clebsch Gordon coefficients (multiplied by 
trivial constants) , and the radial integrals are performed 
numerically. The number of basis functions needed to 
converge the ground state energy to a given relative ac- 
curacy strongly depends on the interaction strength con- 
sidered. At unitarity, e.g., we need a larger basis set 
than in the regime where a s is small and positive (see 
Sees. MM and MM- 

The computational effort of diagonalization schemes 
such as that outlined above increases dramatically with 
increasing number of particles, and eventually becomes 
computationally unfeasible. For larger number of parti- 
cles, we thus resort to an alternative numerical approach, 
the FN-DMC method [5l|, l52j , which exhibits a more fa- 
vorable scaling with increasing number of particles. Our 
implementation of the FN-DMC method has been dis- 
cussed in detail in two recent papers 0, |46[; here, we 
only review the key points. 

The FN-DMC technique, as used throughout this pa- 
per, determines an approximate energy of the many-body 
system whose corresponding eigenfunction has the same 
symmetry as a so-called guiding function ipx, i.e., the FN- 
DMC technique determines the energy of a state that has 
the same nodal surface as ipT but that may differ from 
ipT in other regions of the configuration space. If the 
nodal surface of ipT coincides with that of the true eigen- 
function, then the FN-DMC method results — within the 
statistical uncertainty that stems from the stochastic na- 
ture of the approach — in the exact eigenenergy. If the 
nodal surface of 4>t differs from that of the true eigen- 
function, then the FN-DMC method results in an upper 
bound to the true eigenenergy whose eigenstate has the 
same symmetry as ipx- For example, ipx can be con- 
structed so as to obtain an upper bound for the lowest 
eigen energy with total angular momentum L to t = or 
1 [401 . In this paper, we restrict our FN-DMC calcula- 
tions to the energetically lowest-lying gas-like state of the 
system. 

We consider three different parametrizations of the 
guiding function ipx'- (i) A guiding function ipxi whose 
nodal surface is constructed by anti-symmetrizing a pair 
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function. If N is odd, a single-particle orbital is added 
(with the proper anti-symmetrization) . The detailed 
functional form of 0ti is given by Eqs. (35)-(39) of 
Ref. (46j . (ii) A guiding function ipx2 whose nodal sur- 
face coincides with that of the non-interacting ideal-gas 
nodal surface for the same number of fermions of species 
1 and species 2. The parametrization follows that given 
by Eq. (40) of Ref. j46|. (iii) A guiding function t/> T3 
whose functional form allows, at least in principle, to 
interpolate between the nodal surfaces of -0ti and ipT-2- 
The functional form is given by Eqs. (3)- (4) of Ref. [43| . 

III. RESULTS 
A. Small negative s-wave scattering length 



TABLE I: Energies E NltN2 of the non-interacting system and 
dimensionless coefficients Cn x ,n 2 that determine the ground 
state energy of weakly-attractive trap- and mass-imbalanced 
two-component Fermi gases for selected Ni and N2 values. 
The subscript pair (j, k) can take the values (1, 2) or (2, 1). 
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This section considers the properties of small two- 
component Fermi systems with unequal masses and 
unequal trapping frequencies in the weakly-attractive 
regime, where \a s \ is small (a s < 0). In this regime, a 
compact expression for the ground state energy of the 
Hamiltonian given in Eq. |T]) can be determined within 
first order degenerate perturbation theory for the Fermi 
pseudo-potential Vs , Eq. ^ . Denoting the energy of the 
non-interacting system with Ai atoms of mass m\ and A2 
atoms of mass mi by E^l N2 (see Table U for selected val- 
ues), the perturbative expression for the energy Ejsr lt N 2 
reads 



E 



Ni,N 2 



E 



NI 

N U N 2 



huj- — CnlN 2 , 
a-ho 



where 



and 



miLUl + m2L02 



ah, 




(14) 



(15) 



(16) 



The quantities Q and a^ have been defined so that the 
coefficient C± t i is constant (i.e., independent of 77, see 
below). We refer to uj and (iho as the "natural angular 
trapping frequency" and the "natural oscillator length" 
of the two-body system in the BCS regime. 

The coefficients Cm 1 ,n 2 are listed in Table U for se- 
lected Ni and N2 combinations with \Ni — A2I = or 
1 (N < 8). They reduce to those reported in Ref. [46| 
for equal masses and equal frequencies. The coefficients 
Cni,n 2 i n Tableware written in terms of aho,i and a,ho,2] 
alternatively, they can be written in terms of 77, 



77 = 1 — 



aho,2 
ahoA 



(17) 



The quantity 77 measures the density imbalance of the 
non-interacting two-component Fermi gas. For 7/ = 0, 



z 2 
z 
U 

< 1 
u 1 





1 1 1 1 1 1 1 








* N, 















1 


-0.5 





0.5 ] 



FIG. 1: (Color online) Dimensionless coefficients Ci,i (solid 
line), C2,i (dashed line), Ci,2 (dotted line) and 6*2,2 (dash- 
dotted line) as a function of r\ for weakly-attractive two- 
component Fermi gases. The dash-dash-dotted line shows 
the dimensionless quantity Ca, Eq. (|20[) , which determines 
the excitation gap A(A) for N = 3, 5 and 7. 



the oscillator lengths aho,i [i = 1 and 2) coincide; for 
closed shell systems with Ai = N2, this implies fully 
overlapping densities of the two non-interacting compo- 
nents. For 77 < 0, we have a,ho,2 > a/10,1, while for 77 > 0, 
we have a,ho,2 < dho,i- For k = 4, e.g., 77 < corresponds 
to 0J2 < uji/A and 77 > corresponds to 0J2 > ^i/4. 

Figure Q] shows the dimensionless coefficients Cm 1 .n 2 
for A < 8 as a function of ?/. Plotted this way, the co- 
efficients Cn 1 ,n 2 for fixed N-y, N2 and 77 but different k 
collapse to a single curve. For A — 2, Cn ± ,n 2 is constant 
(see above). For N = 4, Cn ± ,n 2 is maximal for 77 = and 
decreases as |t7| increases. This implies that the attrac- 
tive interspecies scattering length a s can most effectively 
introduce correlations that lead to a lowering of the en- 
ergy, compared to N , when the densities of the two 
components overlap fully. For A = 3, the coefficient C2,i 
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decreases with increasing 77 while the coefficient in- 
creases with increasing rj. This "asymmetry" can be un- 
derstood by realizing that the maximal density overlap of 
the two components for odd-iV systems occurs for finite 77 
and not for -q — 0. For N\ = 2 and A^ = 1, e.g., the maxi- 
mal density overlap of the non-interacting system occurs 
for -q < 0; consequently, the C2.1 coefficient decreases 
with increasing rj. For N± = 1 and N2 — 2, in contrast, 
the maximal density overlap of the non-interacting sys- 
tem occurs for rj > 0. This explains the reversed behavior 
of C'2.1 and Ci 2 as a function of rj. 

The energies for systems with even and odd total num- 
ber of atoms determine the excitation gap A(N) (see, 
e.g., Ref. 



A(N) 



E(N-l)/2,(N+l)/2 + £'(Ar+l)/2 : (A r -l)/2 

2 

^(iV-l)/2,(JV-l)/2 + -E'(A'+l)/2,(A r + l)/2 



> (18) 



where we have taken iV to be odd. Using the perturbative 
energy expression, Eq. (fT4|) . we find 



A(N) 



-hjj- — Ca, 

a-ho 



where 



1 



5a L,i< 



'ho, 2 



4V27T 



(19) 



(20) 



'■ho 



for A" = 3, 5 and 7. The excitation gap determined per- 
turbatively is independent of N for N < 7 for all 
and Wj combinations. A dash-dash-dotted line in Fig. [1] 
shows Ca as a function of 77. The coefficient Ca, and con- 
sequently also A (AT), is largest for rj = and decreases 
with increasing |^|. This can be readily understood by 
realizing that the energies for odd- A" systems [first term 
on the right hand side of Eq. ([T5)l ] average to a constant, 
and that the average of the energies for even-Af system 
[second term on the right hand side of Eq. (fT8"|) ] is mini- 
mal for rj = 0. 

The equal-frequency systems with mass ratio K corre- 
spond to rj — (k — 1)/k. Figure Q] shows that the co- 
efficient Ca decreases with increasing mass ratio k for 
systems with u>i = 0J2, in agreement with the findings of 
Ref. \M- 



B. Infinitely large s-wave scattering length 

This section considers infinitely strongly interacting 
two-component Fermi systems with diverging s-wave 
scattering length a s and varying mass and frequency ra- 
tios. Throughout this section, we express energies in 
units of the average oscillator energy fuv, 



hw\ + hu>2 



(21) 



this unit is convenient since the energies of the non- 
interacting systems with Aq = N2 are directly pro- 
portional to fku (see, e.g., Table U for small N). For 
A" = 3 — 14 atoms, we determine the eigenenergies of 
the stationary Schrodinger equation by the FN-DMC 
method. For A" = 2, we compare the diffusion Monte 
Carlo (DMC) energies (in this case, the ground state of 
the system is nodeless and no nodal approximation needs 
to be made) with the energies obtained from the diago- 
nalization scheme. 

Table HI] reports selected two-body energies (Aq = 
N2 = 1) for k = 1, 4 and 8 at unitarity. The ener- 
gies in the third and fourth column are calculated for the 
square well potential with Rq = 0.01a^ o .i using the DMC 
and the diagonalization approaches, respectively. We an- 
alyzed the convergence of the energies obtained by the 
diagonalization approach by considering basis sets with 
up to about 1000 basis functions. Within the statistical 
uncertainties of the DMC energies, the values reported in 
columns three and four agree. To estimate the energy's 
dependence on the range Ro of the square well potential, 
we diagonalize the Hamiltonian matrix for different Rq. 
We find that the energies for fixed ratio and frequency 
ratios vary linearly with Rq, allowing for a simple linear 
extrapolation to the Ro — > limit (see Ref. [46| for a sim- 
ilar analysis of equal-frequency systems). The energies 
for Rq = 0.01a^ Oi i are slightly larger than the extrap- 
olated zero-range energies (seventh column of Table |TTJ) 
for all mass and frequency ratios considered, and deviate 
by less than 0.5% from the extrapolated zero-range ener- 
gies. Our extrapolated two-body energies at unitarity for 
equal frequencies equal 2ftw for all k, in agreement with 
analytical results for the zero-range potential 57]. For 
larger systems (see below) we do not explicitly extrap- 
olate to the zero-range limit. Based on our two-body 
results, we estimate that the finite range effects of the 
FN-DMC energies for the larger systems at unitarity are 
at most about a few times larger than the statistical un- 
certainties. 

Figure [^a) shows the two-body energies Ei t i calcu- 
lated by the DMC method for the square well potential 
with Rq = 0.01a^ Oj i and l/|a s | = for k = 1,2,4,6 and 
8 as a function of 77. For equal masses, Fig. [2ja) shows 
the energies for frequency ratios ranging from 1/2 

to 1. For unequal masses, the ratio o^/wi of trapping 
frequencies shown ranges from values a bit smaller than 
to 1. In units of fko, the two-body energies for a 
fixed 77 decrease with increasing mass ratio k. Further- 
more, the minimum of the E± t i curves moves to larger 77 
as k increases. 

To shed further light on the behavior of the two-body 
energies, the fifth and sixth columns of Table HT1 show the 
expectation value of H cm , i.e., E^l with (NL) — (00), 
and the ground state expectation value of H re i for the 
square well potential with Ro = 0.01ah o ,i- The sum of 
these two expectation values coincides with the energy 
obtained for a single basis function [namely, <&NLM L 4>nimi 
with (NLMi,nlmi) = (000000)] in the diagonalization 
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TABLE II: Selected expectation values, in units of hui, for 
the two-body system in the ground state with Ni — N2 = 1 
at unitarity for k = 1, 4 and 8 for various frequency ratios 
The energies in column 3 [superscript (1)] are calcu- 
lated for the square well potential with Ro = 0.01ah Oj i using 
the DMC method; in this case, the statistical uncertainty is 
in the last digit reported (or smaller) . The expectation values 
in columns 4-6 [superscript (2)] are calculated for the square 
well potential with Ro = 0.01ah Oi i using the diagonalization 
scheme. The energies in column 7 [superscript (3)] are ob- 
tained by extrapolating the energies obtained by the diago- 
nalization scheme for various Ro to the Tio — > limit; the 
extrapolation error is estimated to be at most O.OOlhcu. 
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1.859 




1/4 


1.867 


1.866 


1.200 


0.725 


1.863 




3/20 


1.918 


1.921 


1.218 


0.784 


1.918 


8 


1 


2.003 


2.003 


1.500 


0.503 


2.000 




3/4 


1.898 


1.898 


1.340 


0.560 


1.895 




1/2 


1.782 


1.783 


1.155 


0.642 


1.780 




1/4 


1.700 


1.701 


0.980 


0.761 


1.697 




1/8 


1.726 


1.724 


0.943 


0.843 


1.720 



approach. The difference between the fully converged en- 
ergies (column 4 of Table [TTJ) and this sum is due to the 
coupling between the center-of-mass and relative degrees 
of freedom. The expectation value of V coup vanishes or is 
negative for all two-body systems considered in this work 
and its magnitude increases for a fixed k with increasing 
lui — LU2- For k = 1, the increase of (H cm + H re i)/(Hu>) 
with increasing lj\ — L02 is larger than the decrease of 
(V coup ) / (ftw)] consequently, the equal- frequency system 
has the lowest energy. For k — 4 and 8, (H cm ) / (fuv) first 
decreases with increasing lu\ — u>2 and then increases for 
UJ2/UJ1 < 1/k, while the quantity (H re i) / (hw) increases 
with increasing lo\ — L02 for all 102/1^1- It can be deter- 
mined readily that the energy of the two-body system 
at unitarity in the zero-range limit without the coupling, 
(uj re i/2 + 3oj cm /2)/w, is minimal at 77 ~ 0.43 and 0.54 
for k = 4 and 8, respectively. Since the absolute value of 
(Vcoup) is fairly small compared to that of (H cm + H re i), 
the minimum of the energy Ei^i/(huj) shifts only slightly 
when the coupling term V CO up is included [see Fig.^a)]. 

The two-body system with unequal frequencies has 
been discussed previously by a number of groups. The 
energies of the lowest-lying gas-like states and the most 
weakly-bound molecular states of the trapped 40 K- 87 Rb 
dimer have, e.g., been measured experimentally and been 
determined theoretically as part of a project on Fermi- 



Bose mixtures in a lattice [53j. Also, the effect of the 
coupling between the center-of-mass and relative mo- 
tions has been investigated in the context of confinement- 
induced resonances [H, [5!^ [6(| • Our main focus lies in 
extending the two-body study presented above to larger 
unequal-frequency systems with three, four or more par- 
ticles per lattice site. While an increasing body of litera- 
ture exists for larger equal-frequency systems, the regime 
where the center-of-mass motion does not decouple has, 
to the best of our knowledge, received only little atten- 
tion for larger systems, despite its immediate relevance 
to ongoing experiments. 

Figures dfb) and (c) show the FN-DMC energies, in 
units of hio, for two-component Fermi gases with N — 
3 and 4 as a function of 77 for various mass ratios n. 
The overall behavior of the N — 4 energies is similar to 
that of the N = 2 energies. For a given 77, the N = 4 
energies decrease with increasing k. For k = 1, the energy 
i?2,2 is minimal for 77 and increases with increasing 
\f]\. As K increases, £"2,2 is minimal for positive 77 values. 
As discussed above, the N = 2 eigenenergies for equal 
frequencies at unitarity approach 2ftw in the zero-range 
limit for all k. For N = 4, in contrast, the energies 
at unitarity depend on the mass ratio even when the 
trapping frequencies coincide [4q |. 

The N = 3 energies at unitarity behave qualitatively 
different from the N = 2 and 4 energies. For systems 
with one spare heavy atom [symbols connected by dot- 
ted lines in Fig. Hfb)] , the energy for a given k decreases 
with decreasing 77. For systems with one spare light atom, 
the behavior is reversed, i.e., the energy increases with 
decreasing 77. This behavior is similar to that of the coef- 
ficients C\$ and 62,1 (see Fig. [T]), and can qualitatively, 
as in the perturbative regime, be explained in terms of 
the overlap of the densities of the two components. 

Table IIIII summarizes selected FN-DMC energies for 
N = 3 — 6 at unitarity calculated for the square well 
potential with Rq = O.Olaho.i- The N = 4 energies for 
77 = and k = 4 and 8 are slightly lower than those 
reported in Ref . [45[ . This is due to the fact that the k > 
1 calculations of Ref. [45| for even N were restricted to 
the guiding function iprii while this paper also considers 
the guiding function i/)t2, whose nodal surface coincides 
with that of the non-interacting system. For N = 4 and 
77 = 0, the energies obtained using ipT2 are lower than 
those obtained using ipTi ■ Table IIIII also indicates the 
total angular momentum L to t of the lowest energy state 
obtained by the FN-DMC method. For N = 4 and 6 (no 
superscript in Table ITTT]) . we obtain the lowest energy for 
a guiding function with L to t — for all mass and trapping 
frequency ratios considered. For N — 3, in contrast, 
the total angular momentum of the lowest energy state 
depends on the system considered and is either or 1. 
For N — 5, the lowest energy state found by the FN- 
DMC method has L tot = 1. 

To investigate the "angular momentum crossover" of 
the three-particle system with two light atoms and one 
heavy atom in more detail, Fig. [3] shows the FN-DMC 
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FIG. 3: (Color online) Ground state energies for N = 3 (two 
light atoms and one heavy atom) in units of huj as a function 
of ft at unitarity; the energies are calculated by the FN-DMC 
method for the square well potential with Ro = 0.01aho,i- 
The symbols connected by dotted and dashed lines are cal- 
culated for equal trapping frequencies and equal trapping 
lengths, respectively. The energies shown by circles and 
crosses correspond to states with total angular momentum 
Ltot = 1 and 0, respectively. 




FIG. 2: (Color online) Ground state energies in units of hio 
for two-component Fermi gases at unitarity as a function of 
r\ for (a) TV = 2, (b) TV = 3 and (c) TV = 4 for various mass 
ratios k for the square well potential with Ro = 0.01a,h o ,i- The 
energies for N — 2 are calculated by the DMC method and 
those for TV = 3 and 4 by the FN-DMC method. The energies 
for k = 1, 2, 4, 6 and 8 are shown by circles, crosses, squares, 
diamonds and triangles, respectively. For a given k, symbols 
are connected by lines to guide the eye: Dotted lines are used 
for TV = 2, TV = 3 with two heavy atoms and one light atom, 
and TV = 4. Solid lines are used for TV = 3 with two light 
atoms and one heavy atom. For k = 1 and TV = 3, the dotted 
(solid) line connects the energies for systems in which two 
particles feel the smaller (larger) trapping frequency. For each 
k, the energies for the largest r\ value considered correspond 

tO LOl — L02- 



energies for equal frequencies (symbols connected by a 
dotted line) and for equal trapping lengths (symbols con- 
nected by a dashed line) as a function of the mass ratio 
k. Energies of states with L to t = 1 and are shown 
by circles and crosses, respectively. For equal trapping 
frequencies, the angular momentum of the lowest energy 
state is L to t = 1 for all mass ratios considered. For equal 
trapping lengths, in contrast, L to t changes from 1 to 
at k i=a 4. This angular momentum crossover has also 
been observed in calculations that employ the correlated 
Gaussian (CG) approach [611 ]. 

Figure 0] shows the energies for TV = 2 — 6 with n = 4 
and 8 as a function of 77. The energies for TV = 2, 4 and 
6 behave similarly as a function of 77, with the TV = 2 en- 
ergies being nearly constant [see also Fig. Ufa)] and the 
TV = 4 and 6 energies showing a stronger decrease than 
the TV = 2 energies as 77 decreases from 3/4 to sa 0.4 for 
k — 4 and from 7/8 to w 0.5 for k — 8 [see also the 
discussion in the context of Figs.[DJb) and (c)]. Further- 
more, the energies E 2y i and £3,2, and the energies £1,2 
and E2.3 show a similar overall behavior. The ordering 
of the energy levels for equal trapping frequencies (for 
k = 4, e.g., this corresponds to 77 = 3/4) is, from bottom 
to top, Ei,!, E ly2 , E 2 ,i, E 2 ,2, E 2 ,3, £3,2 and £3,3. This 
ordering changes as i] decreases; for k = 4, e.g., £2,1 be- 
comes larger than £2,2 at rj « 0.6, and also larger than 
£2,3 at 77 sa —0.2. Similarly, £3.2 becomes larger than 
£3,3 at 77 sa 0.4. 

The small TV energies can be combined to predict the 
energetically most favorable configuration of an optical 



TABLE III: Selected FN-DMC energies E Ni ,n 2 , in units of 
hu, for N = 3- 6 (| JVi - 7V 2 1 < 1) at unitarity for k = 4 and 8. 
The FN-DMC energies are uncertain in the last digit reported. 
The guiding functions used to obtain the energies marked by 
a superscript "*" have total angular momentum Ltot ~ 1 
and those not marked by a superscript have L to t = 0. The 
superscript "?" marks systems for which guiding functions 
with Ltot = and 1 result in energies that are indistinguish- 
able within the statistical uncertainties. For comparison, the 
CG approach results in £2,1 = 5.67/ia), E1.2 = 1.96hu> and 
£2,2 = 4A5huj for k = 8, 012/^1 = 1/8 and a range compara- 
ble to that employed in the FN-DMC calculations [61]; as in 
the equal- frequency case [U 0, [4(|, the FN-DMC energies 
compare favorably with the energies calculated by the CG 
approach. 
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6.88* 
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8.81* 


5.72* 


8.24 


3/20 
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4.50* 
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8.25* 
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2.72* 


4.81 


8.31* 


5.90* 
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5.44 


2.18* 


4.46 


8.56* 


5.08* 


7.90 


1/8 


5.80 


1.96* 


4.46 


8.96* 


4.83* 


7.82 



lattice with small tunneling amplitude, approximately 
harmonic lattice sites, twice as many particles of one 
species than of the other and a filling factor equal to 
or smaller than 3/2. For equal masses and equal trap- 
pin g fr equencies, it has been shown previously for all 
a s |4ll . ml [62| that it is energetically more favorable 
for one spin-up and one spin-down atom to occupy one 
lattice site and for the second spin-down atom to oc- 
cupy a different site [we refer to this as the "(2+1)- 
configuration"] than for the one spin- up and the two spin- 
down fermions to occupy the same lattice site [we refer 
to this as the "(3+0)-configuration"]. Figure O extends 
this analysis to trap- and mass-imbalanced systems at 
unitarity. A dotted line shows the energy of the (2+1)- 
configuration for k — 1 at unitarity as a function of rj [the 
particle in the singly-occupied lattice site feels the larger 
(smaller) frequency for 77 > (77 < 0)], while squares 
show the energy of the (3+0)-configuration. No cross- 
ing between these two curves is observed, with the en- 
ergy of the (3+0)-configuration being larger than the en- 
ergy of the (2+l)-configuration. We find a similar be- 
havior for k = 4: For the frequency ratios considered, 
the energy increase due to placing two like fermions in 
the same lattice site is so large that it is more favorable 
to place the like fcrmion in a separate lattice site instead 
(and thereby "loosing" the energy decrease due to hav- 
ing two like fermions interact with the unlike fermion 
through an attractive two-body potential). For k = 8, 
however, the behavior is different: Figure O shows that 
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FIG. 4: (Color online) Ground state energies at unitarity 
for (a) K = 4 and (b) k — 8 as a funcion of rj for N = 2 
(circles), N = 3 (squares and crosses), N — 4 (diamonds), 
N — 5 (asterisks and triangles) and N — 6 (pluses). Energies 
for even N systems are connected by solid lines, and those for 
odd N systems with a spare heavy and a spare light particle 
by dashed and dotted lines, respectively. The calculations 
are performed using the DMC (N = 2) and the FN-DMC 
(N = 3 — 6) method for the square well potential with Ro = 
0.01a ho 1 . 



the (2+l)-configuration (solid line) is energetically favor- 
able for larger 77 and the (3+0)-configuration (circles) is 
energetically favorable for smaller 77. The crossover is pre- 
dicted to occur at 77 « 0.8. This suggests that it might be 
possible to introduce a macroscopic phase transition of an 
optical lattice system by changing the trapping frequency 
felt by one of the species if the mass ratio is sufficiently 
large but not necessarily so large that three-body bound 
states with molecular character exist. 

We now discuss the odd-even oscillations of larger 
trapped unequal-mass systems at unitarity. In partic- 
ular, we focus on systems with 77 = 0, i.e., on systems for 
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FIG. 5: (Color online) Comparison of the energies at uni- 
tarity for the (3+0)-configuration calculated by the FN-DMC 
method (symbols) and those for the (2+l)-configuration cal- 
culated by the diagonalization scheme (lines) as a function 
of r\. The energies for the (3+0)-configuration are shown by 
squares, asterisks and circles for k = 1, 4 and 8, respectively, 
while those for the (2+l)-configuration are shown by a dot- 
ted, dashed and solid line for n = 1, 4 and 8, respectively. For 
k = 1, two particles feel the larger (smaller) angular trapping 
frequency for r\ > (r\ < 0); for k > 1, the systems considered 
consist of two heavy fermions and one light fermion. 



1.25 



a, 

H 
W 

S w 0.75 

w 



0.25 



1 1 1 1 1 1 

- § 


1 i 1 i 1 i 


e 

- + + ? 


9 9 A3 - 


* 




¥ * * 


(b) - 




8 
N 



10 12 14 



which the two oscillator lengths aho,i and a/ lo ,2 coincide. 
Equal-mass systems with both even and odd N have al- 
ready been discussed in Refs. B, [46| and unequal-mass 
systems with even N in Ref. [45(]. Results for unequal- 
mass systems with odd N (N > 5), in contrast, have not 
been presented before. For even N, N > 6, we obtain the 
lowest FN-DMC energy for K = 1, 4 and 8 for the guid- 
ing function ipTi (see Sec. Ill Bp . whose nodal surface is 
constructed by anti-symmetrizing a two-body pair func- 
tion. For odd N, we find that the guiding function that 
results in the lowest FN-DMC energy depends not only 
on N but also on n: For k = 1, the guiding function ipT2 
results in the lowest FN-DMC energy for N < 9, i^tz 
for N = 11 and ip T i for N > 13. For k = 4 and 8 sys- 
tems with a spare light atom, the guiding functions ipT2 
and ipTi result in the lowest FN-DMC energy for N < 7 
and N > 9, respectively. For n — 4 and 8 systems with 
a spare heavy atom, in contrast, the guiding functions 
■0T2 and ipT3 result in the lowest FN-DMC energy for 
N < 9 and N > 11, respectively. The density profiles 
reveal that the spare particle of the odd-iV systems with 
N > 11 and k = 1 is located predominantly near the edge 
of the cloud 0, [4(|. Thus, one may consider the core 
region as "fully paired" and the edge region as "partially 
paired" . For k > 1 systems with a spare heavy fermion, 
the pairing function ipxi results in a higher energy than 
ipT3 for N > 9, suggesting that the nodal surface ipT3 
allows for a higher probability of three fermions to be 



FIG. 6: (Color online) Energies in units of (a) hui and (b) 
gNi.ETF a ^ unitarity for r\ — as a function of N for k = 
1, 4 and 8. The energies are calculated for the square well 
potential with Ro = 0.01ah o ,i using the DMC (N = 2) and 
FN-DMC (N = 3 - 14) methods, (a) k = 1 (pluses connected 
by a solid line), k = 4 [crosses for even N, and diamonds 
(spare heavy atom) and squares (spare light atom) for odd 
N], and k = 8 [triangles for even N, and asterisks (spare 
heavy atom) and circles (spare light atom) for odd N]. (b) 
K = 1 (uppermost solid line for even TV and pluses for odd 
N), n — 4 [middle solid line for even N, and diamonds (spare 
heavy atom) and squares (spare light atom) for odd N], and 
k — 8 [lowermost solid line for even N, and asterisks (spare 
heavy atom) and circles (spare light atom) for odd N]. 



in close proximity than the nodal surface of ipxi (recall, 
for sufficiently large k three-body bound states with one 
quantum of angular momentum exist) . 

Figures[S{a) and (b) show the energies for systems with 
up to N — 14 atoms with k = 1, 4 and 8. The even-TV 
energies, shown in units of fku in Fig. EJa), decrease for 
a given N with increasing k [for N = 2 and 4, see also 
Figs.[2]Ja) and (c)]. The energies for systems with a spare 
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heavy particle (diamonds and asterisks for k — 4 and 8, 
respectively) are notably smaller than the corresponding 
odd TV energies for k = 1. The energies for systems 
with a spare light particle (squares and circles for k — 
4 and 8, respectively), in contrast, are higher than the 
corresponding energies for k = 1 for small TV and nearly 
coincide with the corresponding energies for k = 1 for 
larger TV. Further optimization of the nodal surface of the 
guiding functions employed in the FN-DMC calculations 
may result in tighter upper bounds for the energies; a 
more detailed investigation of larger odd- TV systems with 
re > 1 is relegated to the future. 

To see the odd-even oscillations more clearly, Fig. |6]Jb) 
scales the energies from panel (a) by the "smoothed" 
extended Thomas-Fermi energies E NI ' ETF of the non- 
interacting system [gU , 



E 



NI,ETF 



= hu 



(3TV) 4 / 3 



1 



(3TV) 



-2/3 



(22) 



For k = 1, the scaled energies follow two distinct curves; 
the curve for odd TV (pluses) is higher than that for even 
TV (topmost solid line) , reflecting the non- vanishing exci- 
tation gap at unitarity (see below and Refs. [43l. H4. |46| ) . 
For fixed TV, the scaled even-TV energies decrease with 
increasing k (45J. The scaled energies for odd- TV sys- 
tems with one spare heavy atom are lower than the corre- 
sponding scaled energies for systems with TV — 1 fermions 
and the same re. The scaled energies for odd- TV systems 
with one spare light particle, in contrast, are higher than 
the corresponding scaled energies for systems with TV — 1 
fermions and the same re. 

Next, we combine our energies for even and odd TV to 
determine the excitation gap A(TV), Eq. ([!§]) , at unitar- 
ity. Figure [7] shows A (TV) for TV = 3 as a function of 
r\ for various re. Although the FN-DMC energies them- 
selves provide an upper bound to the exact eigenenergies, 
the excitation gap is not variational. Figure [7] shows that 
the excitation gap A(3) for re — 1 through 8 collapse — for 
the systems considered — approximately to a single curve 
for j] > 0; A(3) is maximal around r\ = and decreases 
with increasing rj. While more detailed calculations may 
reveal the dependence of A (3) on re for a given rj in more 
detail, our calculations suggest that A (3) is determined 
predominantly by rj and only secondarily by re. This is 
in contrast to the energies themselves (see, e.g., Fig. ^ 
and can be attributed at least partially to the fact that 
the energies of the two odd-TV systems (that with a spare 
heavy and that with a spare light particle) enter into 
A(3) as an averaged quantity. 

For larger TV, we determine the excitation gap A(TV) 
at unitarity for k = 1, 4 and 8 for systems with equal 
trapping lengths, i.e., for ah ,i = aho,2- Figure [5] shows 
that A (TV) , expressed in units of fko, is nearly constant (« 
Q.7hhw) for TV = 3, 5 and 7 for all mass ratios considered 
(for TV = 3, see also Fig. [7]). For TV > 9, the excitation 
gap A (TV) is largest for re = 1 and smallest for re = 8. 
Despite the fairly large uncertainties of the excitation 
gap (see caption of Fig. [5]), we are quite confident that 
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FIG. 7: (Color online) Excitation gap A(3) in units of hu 
as a function of r\ for two-component Fermi gases at unitarity 
for various re. The A(3), calculated using the DMC (TV = 2) 
and FN-DMC (TV = 3 and 4) energies, are shown by circles, 
crosses, squares, diamonds and triangles for re = 1, 2, 4, 6 and 
8, respectively. Dotted lines connect symbols for a fixed re to 
guide the eye. 
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FIG. 8: (Color online) Excitation gap A(TV) in units of hco 
at unitarity for 77 = (i.e., for miUi = 7712^2) as a function 
of TV for re = 1 (circles with errorbars), re = 4 (squares), and 
re = 8 (diamonds). The errorbars of A(TV) for re = 4 and 8, 
which are not shown to enhance the clarity of the figure, are 
a bit larger than those for re = 1. The excitation gap A (TV) 
is calculated from the energies shown in Fig. [6] 



the excitation gap does indeed decrease with increasing re 
but fixed TV (TV > 9). This is a direct consequence of the 
decrease of the energies for odd- TV systems with a spare 
heavy particle with increasing re. 
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C. Small positive s-wave scattering length 



This section discusses the behavior of two-component 
Fermi gases with unequal masses and unequal frequencies 
in the BEC regime where the s-wave scattering length a s 
is small and positive. In this regime, the behavior of 
the molecular Bose gas is expected to be gove rned by the 
dimer-dimer scattering length add [64l65l| . For equal- 
frequency systems, the four-fermion spectrum has been 
compared with that of two bosons (4f| [46|, validating 
the "dimer picture" . This section extends the previous 
analysis to two-component Fermi gases with unequal fre- 
quencies. 

We determine the two-body energy Ei.i for unequal 
frequency systems using the diagonalization scheme. For 
a sufficiently small atom-atom scattering length, we find 
that the two-body energy is to a very good approxima- 
tion given by the sum of the expectation values of H cm 
and H re i (the expectation value of the coupling term 
V coup is smaller than 10 _4 fiw for the systems considered 
in Fig. [5]). Thus, an approximate but highly accurate 
expression for the ground state energy of the two-body 
system on the BEC side reads 



El, 



(H re l) 000 + 

cm ■ 



(23) 



Assuming diatomic molecules form, the lowest energy of 
the four-particle system with N± = N2 = 2 and small a s 
can be written as 



3 



2,2 



2(H„ 



rei/000 



— fkd cm 



3„ 



rel,boson ; 



(24) 



where the first term on the right hand side is the internal 
energy of the two bosonic molecules, the second term on 
the right hand side is the center-of-mass energy of the 
two-boson system and the third term on the right hand 
side is the relative energy of the two-boson system. We 
rewrite the latter as 



E r 



= —filOr 

2 



3, 



(25) 



and evaluate the "dimer-dimer interaction shift" Edd in 
first order pert urbation theory (assuming a Fermi contact 
potential) [13, [E3, 



3, 



dd 



add 



(26) 



Here, aho,cm denotes the oscillator length associated with 
Lu C m, i-e., aho.cm — \J fc/ (Muj cm ). It follows that the en- 
ergy difference £2,2 — 2£'i i i should be given by the in- 
teraction shift Edd- The reasoning outlined here for four 
fermions can be extended to larger systems: The trapping 
frequency u re i determines — together with the atom-atom 
scattering length a s — the internal binding energy of the 
molecules while the trapping frequency to cm determines — 
together with the dimer-dimer scattering length add — the 
properties of the composite boson system. 
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FIG. 9: (Color online) Interaction shifts in units of ftw cm 
for four-particle system in the BEC regime. Circles and dia- 
monds show the energy difference £2,2 — 2Ei t i for re — 1 and 
4, respectively; the uncertainty of (£2,2 — 2_Bi > i)/(fioj cm ) is 
about 0.01hu} C m- The energies £2,2 and are calculated 
by the FN-DMC and diagonalization approaches, respec- 
tively, for the square well potential with _Ro = 0.01ah Oi i and 
a s — 0.1a; lo ,i. For comparison, solid and dashed lines show 
the interaction shift Edd determined perturbatively, Eq. (126|l , 
for re = 1 and 4, respectively. The quantity a s /aho,cm in- 
creases with increasing U12/UJ1, which explains the increase of 
the interaction shift with increasing ^2/1^1 • 



Figure [5] shows the difference between the FN-DMC 
energy £^2 an d twice the two-body energy E^i as a 
function of the frequency ratio W2/W1 for k = 1 (circles) 
and k = 4 (diamonds). For comparison, solid and dashed 
lines show the interaction shift Edd, Eq. (|2"6"1) , for n = 1 
(using a dd = 0.608a s JM, EE H, 1(|) and k = 4 (using 
add = 0.77a s [ljj, |45|, |46j ) , respectively. The numerically 
determined energy differences are a bit larger than the 
perturbative result, which is in agreement with the fact 
that the FN-DMC energy En. 2 provides an upper bound 
to the true eigenenergy [52j |. The agreement between 
the interaction shift obtained by solving the full four- 
body Schrodinger equation and by treating the weakly- 
interacting two-boson system perturbatively is similarly 
good for all trapping frequencies considered. This con- 
firms that the relevant "boson frequency" is indeed, as 
has been argued previously by others 0, |47j , given by 



IV. CONCLUSION 

This paper determines the ground state energies of 
two-component Fermi systems under external harmonic 
confinement with unequal masses and unequal frequen- 
cies. We considered the weakly- interacting, small \a s \ 
regime with both positive and negative scattering lengths 
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a s as well as the strongly-interacting unitary regime 
where the s-wave scattering length diverges. In all three 
regimes, we identified convenient energy and length units. 
In the weakly-attractive regime, we treated the atomic 
Fermi gas perturbatively. In the unitary regime, we de- 
termined the eigenenergies numerically. In the weakly- 
repulsive regime, we compared numerical results with 
those obtained by treating the weakly-repulsive molec- 
ular Bose gas perturbatively. 

The calculations at unitarity are performed for a short- 
range potential with small range; the resulting energies 
are estimated to be quite close to the zero-range limit. 
We determine the energies as a function of both the ra- 
tio between the masses of the two species and the ratio 



between the trapping frequencies felt by the two species. 
The small N results presented cover a wide range of mass 
and frequency ratios and can easily be extrapolated to ex- 
perimentally relevant parameter combinations (such as 
6 Li- 40 K mixtures). Our FN-DMC energies provide an 
upper bound to the true ground state energy of mass- 
and trap-imbalanced two-component Fermi systems and 
may be serve as a benchmark for other approaches. 

DB is grateful to J. von Stecher for calculating the 
CG energies reported in the caption of Table IIIII and 
for discussions during the initial stage of this work. DB 
also gratefully acknowledges support by the NSF through 
grant PHY-0555316. 
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